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Abstract. We study the energetic aspects of hybrid pressure-gravity modes pulsations. The 
case of hybrid j3 Cephei-SPB pulsators is considered with special attention. In addition to the 
already known sensitivity of the driving mechanism to the heavy elements mixture (mainly the 
iron abundance), we show that the characteristics of the propagation and evanescent regions 
play also a major role, determining the extension of the stable gap in the frequency domain 
between the unstable low order pressure and high order gravity modes. Finally, we consider the 
case of hybrid S Sct-7 Dor pulsators. 



1. Introduction 

In the HR diagram, many families of pulsators appear as couple. In the hot part of the main 
sequence, we encounter the f3 Cep stars with low radial order modes {y m 3 — 8 d" 1 ) having a 
mixed p-g mode nature (p-mode type in most of the envelope and g-mode type in the region 
just above the convective core). And just at their cool side the Slowly Pulsating B (SPB) stars 
with high radial order gravity modes (y ~ 0.2 — 1 d _1 , n ~ 20 — 50). The same occurs at lower 
temperature: we find the 5 Set stars with low radial order p-modes and just at their right the 
7 Dor stars with high radial order g- modes; and at later evolution stages also we find the couple 
of sdB and Betsy stars. We concentrate here on the case of (3 Cep and SPBs. Their modes 
are driven by a K-mechansism operating in the Iron Opacity Bump (IOB) at logT ~ 5.3. At 
first sight, the different pulsation periods and locations in the HR diagram of these two families 
seem to take their origin in the location of this IOB. As T c ^ decreases, it goes deeper into the 
star, entering progressively the domain of g-modes with larger periods. But recent observations 
complicate this simple scenario: hybrid stars with both types of oscillations at the same time 
are observed. The goal of this paper is to explain their existence and show some aspects of their 
high potential for asteroseismology. 



2. Hybrid (3 Cep - SPB stars 

Recent observations show the existence of hybrid B stars with at the same time oscillations of 
(5 Cep type [y « 4 — 8 d -1 ) and high order g-modes of SPB type (u « 0.2 — 1.2 d _1 ). The 
clearest cases are 19 Mon [2], v Eri [9] and 12 Lac [10]. Another very promising star is HD 180642 
(a primary target of COROT). Several other cases exist such as 7 Peg [3], but they must be 
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Figure 1. Frequencies of unstable I = 2 (left) and £ = 1 (right) modes of (3 Cep-type (blue) 
and SPB-type (red) as a function of T e g , along the evolution sequence of 10 M from ZAMS 
to TAMS with OP opacities. 



considered with caution because of binarity or because of the low S/N and the risks of aliases, 
asking for more precise observations. 

From a theoretical point of view, [15], [13] and [14] showed that with OP opacities ([18], [19]), 
the instability strips of j3 Cep and SPBs intersect, so that theoretical models with both unstable 
p-modes and high order g-modes exist. An illustration is given in Fig. 1 for 10 M models. 

For the structure models of this study, we used the stellar evolution code CLES [17], with by 
default the following physical prescriptions: the new OP opacities [19], the OPAL2001 equation 
of state [16], the mixture of elements by [1], with the enhancement of Ne proposed by [7]. 

3. Driving and damping regions 

To identify the driving and damping regions, we give in Fig. 2 the work integrals for different 
modes of a well chosen 10 Mq model. Left panels correspond to I = 2 modes and right panels 
to I = 1 modes. Regions where W increases (resp. decreases) outwards have a driving (resp. 
damping) effect. We analize first the 1 = 2 modes (left). Unstable modes of j3 Cep type (low 
radial order mixed modes) are given in blue in the top left panel and unstable modes of SPB 
type (high radial order g-modes) are given in red in the bottom left panel. For these two types 
of unstable modes, we see clearly the significant driving around the IOB at logT ~ 5.3. But g- 
modes of intermediate radial order (green curves) are stable because of a significant damping in 
the very deep layers just above the convective core. This damping is a classical radiative damping 
mechanism. The short-wavelength oscillations in the g-mode cavity due to the large values of 
the Brunt- Vaisala frequency lead to very large values of the derivatives of the eigenfunctions 
(d/dr(/(r) sin(J Mr)) ~ kf(r) cos(J kdr) with k ~ \fl\l + l)o" _1 / N/rdr » 0). In particular, 
the temperature gradient variations due to the oscillations are very large in this region; it is 
easily seen that they lead to significant loss of heat during the hot phase (see Eq. 7 of [4] ) , which 
always damps the modes. But why this significant damping occurs for some g-modes and not 
for others ? 

4. Propagation and evanescent regions 

Changes in the respective weights of the driving and damping regions originate in the shape of 
the eigenfunctions. To understand this, classical propagation diagrams are very useful, as shown 




Figure 2. Work integrals f™ dW/ J \dW\ for different £ = 2 modes (left) and £ = 1 modes 
(right) of a 10 M model with T eff = 21245 K, log(L/L ) = 4.06 and OP opacities. 
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Figure 3. Square of the dimensionless Brunt- Vaisaia frequency N 2 B?/(GM) (red), 
Lamb frequency Sf R 3 /(GM) (black) and angular pulsation frequencies uj 2 = a 2 R 3 /(GM) 
(horizontal blue lines), for I = 2 modes (left) and £ = 1 modes (right). 



in Fig. 3 for the 10 M model considered previously. According to the asymptotic theory, the 




Figure 4. Real parts of the transversal displacement £/j (blue) and 

A{r)i h = pi/Vc-^ISf/a 2 - l] 1 ^ \N 2 /a 2 - l)" 1 / 4 £ h (red), for the modes I = 2, g 23 (left) 
and £ = 2, gg (right). 



transversal component of the displacement (dominating for g-modes) behaves as: 
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where the function i 7 is a cosine of constant amplitude in the propagation regions (green in 
Fig. 3), and an exponential decreasing outwards in the evanescent regions (magenta in Fig. 3). 
Fig. 4 (left panel) shows that this formula applies well for the mode £ = 2, g23. The propagation 
region is very large for this mode going up to the IOB, the change of amplitude of ^ is dominated 
there by the factor p -1 / 2 which leads to values much smaller near the center compared to the 
surface. Hence, the deep radiative damping is small compared to the K-driving in the IOB and 
this mode is unstable. The right panel of Fig. 4 shows the same eigenfunctions for the mode 
£ = 2, gg. Now the propagation region is smaller with a larger evanescent region above it (see 
Fig. 3). The corresponding exponential decrease leads to amplitudes much smaller around the 
IOB compared to the deep layers. Hence the deep radiative damping overwhelms the Ac-driving 
and this mode is stable. 



5. Dependence with the spherical degree £ 

In the right panels of Fig. 2 the work integrals for the same model as before are given, but now 
for £ = 1 modes. We see in this case that the (5 Cep-type mixed mode I = 1, gi is predicted 
to be unstable (blue, top right panel), but all £ = 1 high radial order g-modes of SPB type are 
predicted to be stable, contrary to the 1 = 2 case. This can be understood as follows. At high 
frequencies, the eigenfunctions do not depend significantly on I. Hence the range of excited 
modes of (5 Cep type (as a function of the frequency) remains essentially the same for different 
I (blue diamonds in Fig. 1 and blue curves in the top panels of Fig. 2). At lower frequencies, we 
have seen that the respective weight of the damping and driving regions is determined by the size 
of the propagation and evanescent regions. As Sf oc £{£ + 1), g-modes having the same values 
of £(£ + l)/cr 2 have the same propagation regions. This corresponds to a lower frequency for 
£ = 1 compared to £ = 2 modes (compare the left and right panels of Fig. 3). Hence the critical 
frequency at which the upper boundary of the propagation region reachs the IOB is smaller for 
£ = 1 modes. At this smaller frequency, the transition region where the pulsation period is of 
the same order as the thermal relaxation time is deeper in the star, around the hot side of the 



IOB; so the damping effect of this region is larger than the driving of the more superficial layers 
(red curves in the right bottom panel of Fig. 2) and all high order g-modes are stable. 

6. OP versus OPAL 

[13] showed that many more hybrid (5 Cep-SPB are predicted when OP opacity tables ([19]) 
are used compared with OPAL ([12]). This can be seen in Fig. 5 where we compare the range 
of excited modes obtained in the two cases for 10 M models. The comparison of the OP 
and OPAL opacities and their derivative d\n.K/d\nP\ s presented in Fig. 6 shows significant 
differences around the hot wing of the IOB. Fig. 7 compares the work integrals obtained with 
OP and OPAL for a typical (5 Cep-type mode (left) and SPB-type mode (right). The results 
are very close for j3 Cep-type modes because the opacities are essentially the same with OP and 
OPAL in the superficial layers (logT < 5.2) where p-modes have significant amplitudes. But the 
results are different for the high order g-modes: with OP they are excited while most are stable 
with OPAL. To understand this, I recall that two things are required for an efficient driving 
of the g-modes. First, the transition region where the pulsation period is of the same order 
as the thermal relaxation time must coincides with the IOB. And second, this period must be 
large enough so that the corresponding evanescent region is negligible. These two requirements 
are fulfilled with OP; but with OPAL the period at which the transition region and the IOB 
coincide is too small, leading to a too large evanescent region. This explains why the driving of 
high order g-modes (probing these deep layers) is more efficient with OP than with OPAL, and 
more hybrid models are predicted with OP. 
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Figure 5. Frequencies of unstable t = 2 
modes as a function of T e ff , along the evolution 
sequence of 10 M from ZAMS to TAMS, using 
OP (blue) and OPAL (red) opacities. 



Figure 6. Opacities (top) and din K,/dln P\ s 
(bottom) obtained with OP (blue) and OPAL 
(red) opacity tables, 10 M models with T e fj = 
21245 K and log(L/L ) = 4.06. 
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Figure 7. Work integrals for the mode £ = 2, g3 (left) and I = 2, g2o (right) of the same 
models as Fig. 6, using OP (blue) and OPAL (red) opacities. 

7. Hybrid 5 Sct-7 Dor stars 
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Figure 8. Propagation diagram (top) and work 
integrals (bottom) for the modes I = 1, P2, g6 
and g25 of a young 1.55 M Q 5 Sct-7 Dor model. 



Figure 9. Frequencies of unstable I = 1 — 4 
modes of 7 Dor type (red) and 5 Set type 
(blue), for the same model as Fig. 8. 



Lower along the main sequence, hybrid models with both unstable 5 Set-type mixed modes 
and 7 Dor- type g-modes are also predicted ([4], [6]) The observational search of such hybrid 
models was less conclusive than in B stars but some cases were found: the Am star HD 8801 
[11] and the binary HD 209295 [8]. The driving of these stars is very different from the one of 
B stars, being complicated by the time-dependent coherent interaction with convection for cool 
models [4]. However, we point out here that some aspects of the damping mechanisms remain 
very similar. 

In Fig. 8, we give a propagation diagram (top, dimensionless frequencies) and the work 
integrals (bottom) for the modes t = 1, P2, g6 and g25 of a young 1.55 M Q 5 Sct-7 Dor hybrid 



model with T eS = 7330 K, log(L/L ) = 0.746 and log g = 4.295. This model is computed 
with OPAL opacities and Z = 0.02. We used for the non-adiabatic computation the Time- 
Dependent Convection (TDC) treatment of [5]. The mode i = 1, P2 (blue, v = 28.8 c/d) is 
a typical unstable mode of 5 Set type; TDC plays a central role in its driving and damping. 
The mode t = 1, g25 (red, v = 0.98 c/d) is a typical unstable mode of 7 Dor type; the small 
driving at logT ~ 5.3 comes from a k- mechanism in the IOB and the significant driving at 
logT ~ 4.8 is a flux blocking mechanism operating a the base of the convective envelope. The 
mode £ = 1, ge (green, v = 3.66 c/d) is a stable mode which is significantly damped in the deep 
radiative g-mode cavity. We note in the top panel of Fig. 8 that the sizes of the propagation and 
evanescent regions for these modes are similar to the case of B stars considered above. Hence, 
the explanation of the stable gap between the 5 Set and 7 Dor unstable modes is the same as 
before: these g-modes of intermediate radial order have a large evanescent region; hence the 
amplitudes are large in the deep g-mode cavity, and the radiative damping occurring there is 
much larger than the driving of the superficial layers. We note that the model considered here 
is very young, without region of variable molecular weight above the convective core. Hence, no 
bump of the Brunt- Vaisala frequency is present there. 

In Fig. 9, all the frequencies of unstable £ = 1 — 4 modes are given for the same model as 
before. We note that the stable gap between the 5 Set and 7 Dor unstable modes decreases as £ 
increases. This is simply due to the fact that Sf oc £{l + 1), so that the size of the propagation 
cavity of g-modes increases with £ for a given frequency. 

8. Conclusions 

Theoretical models and observations indicate the existence of hybrid models having at the same 
time unstable pressure and gravity modes. Concerning first the driving mechanism of hybrid 
(3 Cep-SPB models we have shown, in addition to its already known sensitivity to the iron 
abundance, that the characteristics of the propagation and evanescent regions also play a major 
role in this context, determining the extension of the stable gap in the frequency domain between 
the unstable low order pressure and high order gravity modes. The same phenomenon occurs for 
the hybrid 5 Sct-7 Dor pulsators. Hence, the comparison between the observed and theoretical 
ranges of unstable modes allows us to constrain the Brunt- Vaisala and Lamb frequencies in the 
very deep interior of these stars. 
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